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Abstract 

We investigate one of the simplest multi-species generalizations of the one dimensional exclu- 
sion process with reflective boundaries. The Markov matrix governing the dynamics of the system 
splits into blocks (sectors) specified by the number of particles of each kind. We find matrices con- 
necting the blocks in a matrix product form. The procedure (generalized matrix ansatz) to verify 
that a matrix intertwines blocks of the Markov matrix was introduced in the periodic boundary 
condition, which starts with a local relation (Arita et al, J, Phys. A 44, 335004 (2011)). The 
solution to this relation for the reflective boundary condition is much simpler than that for the 
periodic boundary condition. 

1 Introduction 

The asymmetric simple exclusion process (ASEP) is a lattice-gas model of interacting particles |16j . 
where each particle is a random walker hopping from a site to one of the neighboring locations only if 
the target site is empty. The ASEP in one dimensional lattice (Z or its subset) has been intensively 
studied as an exact solvable non-equilibrium model [3 [ini Hi] which is related to growth phenomena 
[T51 [TU and applied to modeling of various transport systems "from molecules to vehicles" [T71 [35] . 
One of standard generalizations of the ASEP to multi-species systems (A^-species ASEP) where each 
site i takes a local state ki e {1, . . . , A^-f 1} (A^ > 0) is as follows; nearest neighbor pairs of local states 
(a,/?) = (fci,fci-|_i) are interchanged 

aP^I3a with rate | 1" ^ (1) 
\ q {a> P), ^ ' 

where we impose < g < 1 without loss of generality. We say that the site i is occupied by an 
ath-class particle for ki ^ a < N . We also say that the site i is empty for ki — N + 1, which is 
usually denoted by 0. The usual ASEP corresponds to A^ = 1. We will be formally concerned with the 
zero-species ASEP (A^ = 0) as well. We will call the cases q — 1 and q — multi-species symmetric 
simple exclusion process (SSEP) and multi-species totally ASEP (TASEP), respectively. 

In [5] , the spectral structures of Markov matrices (which govern the dynamics of the system) were 
clarified. On the other hand, in [T^, the matrix product form for the stationary state was found in 
the periodic-boundary case. Then these two studies were combined in 0] where matrices connecting 
dynamics of different values of A^ were constructed. This generalizes the procedure of the matrix 
(product) ansatz for stationary states [71[TT1[T2]. In particular this generalized matrix ansatz enables 
us to transfer information to a system consisting of A^-species particles from simpler systems consisting 
of A'''(< A^)-species particles. The question is whether this generalized matrix ansatz is applicable to 
the reflective-boundary case or not. This will be answered positively, which is the main purpose of 
this paper. 
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Here we define the model on the L-site chain precisely. Let {|1), . . . , + 1)} be the basis of the 
single-site space C^"*"^, and represent a particle configuration fci • • • fc^ as the ket vector \ki ■ ■ ■ ki,) = 
\ki) I/cl) S (C^"*"^)®^. We also use the corresponding bra vectors (fc| and (fci---fci| = 

(fell (g) • • • (g) (fcil with {j\k) — Sjk- In terms of the probability vector 

\P{t))^ P{ki---kL;t)\k,---kL), (2) 

l<fei<Ar+l 

with each coefficient P{ki ■ ■ ■ k^^t) representing the probability of finding the configuration ki ■ ■ ■ k^ 
at time t, our model is governed by the master equation ^\P(t)) = Af'-^'|P(t)). The matrix M'^^-' 
has the form 



V / 2,2+1 

l<i<L-l ' 
N+1 

Mi^c ^ 5] (-e(a-/3)|a/3)(a/3|+e(a-/3)|/3a)(a/?|), (4) 

a, 13=1 

where [M^^ ) acts nontrivially on the ith and (z + l)st components of the tensor product, and 
8 corresponds to the transition rate 

r 1 (a</?), 

e{a-l3) = { q {a>p), (5) 
I {a^p). 

We call the {N + 1)^ x {N + 1)^ matrix A/^^^ total Markov matrix or simply Markov matrix, and the 
{N + 1)^ X [N + 1)^ matrix m'^^ local Markov matrix. The relevant two-dimensional vertex model 
is a special case of Perk-Schultz model [H [31 [T51 [^ . 

We emphasize that we will investigate the multi-species ASEP on the L-site closed segment with 
the reflective boundary condition {free boundary condition). I.e., we do not impose any boundary term 
in the Markov matrix ([3]). The total Markov matrix obviously preserves the number of particles of 
each class and thus it has the following block diagonal structure: 

mW=0M„, M^eEndK™, (C^+i)«i=0v;„. (6) 

m m 

Here each sector V^^ is specified by the number of particles of each class 

Kn= C|fci...fci), (7) 

for m = (mi, . . . , niN+i)- This means that each sector is spanned by the ket vectors corresponding to 
permutations of the sequence ■ ■ N + \ ■ ■ ■ N + \. We also call the label m for the sector 

mi m2 niiv + l 

Vm "sector" . We write the dual space for each sector as 

C= C(fci.--fci|. (8) 

This paper is organized as follows. In section[5]we review the symmetry of the model, and construct 
matrices which connect dynamics of different sectors. We call these connecting matrices "conjugation 
matrix" , and in section[3]we will reconstruct them in a matrix product form. In section|4]we consider a 
relation between Markov matrices of different values of N , by introducing similar conjugation matrices. 
In section[5]we investigate the stationary state and the relaxation to it. We will see that the stationary 
state can be regarded as a product of conjugation matrices. Section [S] is the summary of this paper. 
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2 Symmetry 

The Markov matrix is mapped to a Uq{SU (A^+1)) invariant quantum Hamiltonian H by the similarity 
transformation El El m [501 HI] 

v'' l<fei<JV+l 

H^'''- E (<cO„,,^^["i-- E E (10) 

l<i<L-l i<c,ii<N + i l<a.p<N+l 

Thus M^^) also has a symmetry, which we review in this section. We define Uq{SU{N+l)) generators 
/(") = |n+l)(7i|, e(") = |7i)(7i + l|, fc(«) =^|n)(„| + ^|n + l)(n+l|+ ^ (11) 

l<x<N+l 

and the comultiplications 

A/(") = + Ae(") = e(") (8) (fc("))^i +/d«)e("\ A/c^") = /c^") (12) 

which commutes with the local Hamiltonian: 

[A/'"',<e'] = [Ae("),<i] = [Afc(«),<i] = 0. (13) 
From these local relations, we find global commutation relations 

where 

= /I + fcj" Vi"^ + • • • + A:^"^ • • • 1 (15) 

^(") = el"^(A;^"V • • • (fci"V' + • • • + ei"2i(fci"V + ei"\ (16) 
= fcj") • • • (17) 

and each x^"-* (x = /, e, k) acts nontrivially on site i. From the similarity transformation we have 

[^(") ^ = , m(^)] = , m(^)] = 0, (18) 

where = 5'-iX(")S' {X = F,E,K). By the direct calculation, one can show that the elements 
of and are given as 

(ji...ji|F(")|A:i...fci) = <^ j, = fc,(.^z)), (19) 

[ (otherwise), 

, I n«#{»'>*l'=.'="+i} (3« : - 1 = fci = f^, 

- <^ ^'^ j; = fc,(,^z)), (20) 

[ (otherwise), 
Q = g3(#{^'lii'=«}+#{i'b>'="+i}-i). (21) 

We notice that the matrix F^^'' (resp. i?^")) sends a ket (resp. bra) vector in the sector m = 
(mi, . . . , mN+i) to the sector m(") = (mi, . . . , m„_i, m„-l, m„+i+l, m„+2, • ■ • , mN+i), i.e. i^*^")!/™ C 
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V^(n), V*iE^"'' C y*(„). In other words, these matrices change one particle 71 to rt + 1. From the 
commutation relations (|18p . we have 

i?(")M,„ = M„(„,i?("), M„J(") = -B(")M„,(„, , (22) 

which implies Spec(Mm) C Spec(M„(„)) or Spec(Af,„) ^ Spec(M„(„)) according to diml^i < 
dim or dim^m > diml^cn), respectivelyjll Here Spec(Mm) is the multiset of all eigenvalues 

of sector m, where the multiplicity of each element corresponds to the degree of degeneracy. 
For a sector m = (mi, . . . , totv+i), we define 

m^")" = (mi, . . . , m„_i, m„ - /i, to„+i + ^, m„+2, • ■ • , mAr+i), (23) 

i.e. the sector m^^"^^ is obtained by changing /i particles of nth class to (n + l)st-class particles. The 
matrix (F'"^)'' sends a vector in Kn, to l^(Ti)f , and one can show by induction that each element of 
(F("))'' is calculated as 

(j-i---jL|(i^'"^)^|fci---fcL) = M!g^9^'^^*^''<^l^''="^ (24) 
if 3/ C {1, . . . , i} (#/ = fi) such that - I = k, = n (l e I) and j,, = k,{L ^ I), (25) 

or otherwise. Here [/i]! is the q-factorial ni<p'<;i[/^'] with the g-integer [fi'] = l + q + - ■ • + and 
Q is defined by (PT]) . Similarly, (iJ^"))'' sends a vector in V*^ to V^^(„)m , and each element is calculated 
as 

{k^---kL\{E^^^r\n---jL) = [;,]!QMqE.e.#K>^|fc.'="+l} (26) 
if the condition (|25l) is satisfied, or otherwise. From equation ([22) we have 

(F("))'^Af„ = Af,„(„,. (F("))^ A/,„(^("))'^ = (^("))''M,„(„,. . (27) 
We call this type of relations conjugation relation and the matrix that satisfies it conjugation matrix 

m- 

Comment: Since iJ^^"^ (|TOl) is a symmetric matrix, we have Al'^ — S^MS ^. This means our 
process satisfies the detailed-balance condition, and the stationary-state probability is expressed as 
Po{ki---kL) = {ki ■ ■ ■ kL\S~'^\ki ■ ■ ■ kL) = Ei<i<j<i. sign(/ci-fcj)^ r^j^j^ Stationary state can be 

rewritten in terms of a matrix product form, which we will achieve in another way, i.e. by using the 
generalized matrix ansatz, in section [5l 



3 Matrix product interpretation 



In this section, we write the elements of (F^"))^ and (F^"))^ in a matrix product form. We first define 
matrix-valued matrices fe^^"^ and 6'^"^ of size + 1) x [N + 1) as 



&(") = l\x){x\+D\n){n\+A\n + l){nl 

l<jc< Af + l 

&(") = ^ l\x){x\+D\n){n + l\+A\n + l){n + l\, 



l<x<N + l 



where 



/ 1 \ 
1 



D 







V 



A = 



V 



^ For the periodic-boundary case, this spectral inclusion is not generally satisfied [6]. 



(28) 
(29) 



(30) 
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satisfying the relations DA = qAD. We also use vectors {{w\ and defined as 



{w\ = {lOQ ■■■), \v)) = 



(31) 



which satisfy {{w\A — and D\v)) = \v)). 

The elements of (F^"))^ ^ and ^ can be interpreted as follows: for configurations 

ji • • • JL and fci • • • fci in the sectors m and m'"^^ (j23p . respectively, we have 



■b^"l\v): 



(ji---jL|(F(")nfci...fci) = [/i]!Q^((«;|6;.X---.^.^,^ 
(fci...fci|(i?(")nji...j,.) = [/i]!Q^((«;|Bg^...6£'^.J«)), 



(32) 
(33) 



where we write 6^^' = and 6^,"^ = Equation ([32]) is understood as follows. Let 

Ji • • ■ JL and fci • • • fci be configurations such that 

ji = /ci or ji - 1 = fci = n, (34) 

and set 



= n, fj. = ^ - 1 = fcj = n}, = #{i\jt = h = n}. (35) 

n, ^ 

and &j"''„ = ^, D. Noting X^ie/ < ~ corresponds to how many times we need to exchange 

DA — > AD for reordering the matrix product • • • foj"' „ to A ■ ■ ■ A £) ■ ■ ■ Z3 , we find 



Since = 1 for h ^ we have {{w\h^^ll ■ ■ ■ bf^, H - • • • 



If the configurations do not satisfy the condition ([M)) . we have • • • 



0. Noting this 



and equations (|24l) and p6p we find the matrix product interpretation p2p . One can show equation 
P3p in the same way. Defining the matrices 



we rewrite the conjugation relation (j27p as 



(«) 



(37) 
(38) 

(39) 



In what follows, we show these relations in another way starting from local relations different from 

Using the relation DA = qAD, one can check that the tensor products 6^"^ (8) and 6^"-' (8) 
commute with the local Markov matrix: 



This leads to a global commutation relation 



= 0. 



= 0. 



(40) 



(41) 
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Bookending each element between vectors ((wj and \v)) (PT|) . we obtain (scalar-valued) matrices of size 
(N + l)^ X {N + 1)^ 



\1>(' 



(42) 



which satisfy 



= 0. 



(43) 



Since ipm,m("')^' f37)) and ifmi")'' ,m (EH) are submatrices of ^^"^ and , respectively, the commutation 
relations (|43l) lead to the conjugation relations ((39|) . 



4 Relation between dynamics of different values of N 

In this section, we construct a matrix which intertwines dynamics of an A'^-species sector m = 
(mi, • • • , nipf+i) and an {N — l)-species sector m' = (mi, . . . , mn-i,rnn + m„+i, m„_|_2, • • • , tow+i)- 
The dynamics of the sector m' is essentially same as that of the sector m = (mi, . . . ,m„_i,m„ + 
m„+i, 0, m„+2, ■ • • , TTT-N+i) by regarding particles of xth class (x > n + 2) as that of (x— l)st class. 
(Note that 

= jji '^ Thus we have already known that there exit matrices 4'mm' and (pm'm 
that satisfy conjugation relations Mm^pmm' = ^mm>Mm> and ipm'mMm = Mm'ipm'm- 

Now we restrict our consideration to sectors (mi, . . . , mAr+i) such that mi > for all i {basic 
sector), and introduce alternative labeling for basic sectors [B]: 

m = (mi, . . . ,mAr+i) o s = {si, . . . , sat} (1 < si < • • • < sat < L - 1) (44) 

with the correspondence Si — mi + • • • + m^, or equivalently rrii = Si — Si-i (sq = 0, sn+i — L). In 
particular, for the zero-species sector, we use the labeling 

[L) o 0. (45) 

According to this correspondence, the sector m' is labeled by s \ {s„}. We write = Vm and 
— Mm for s f^- m. 

First we construct V'mm' = V's,s\{s„} that satisfies the conjugation relation by starting with the 
following matrix- valued matrix a^^'^'^ of size (A^-|- 1) x A^, a "degenerated version" of the matrix 6^"^ 
(EHl): 



= ^ t\x){xn{x)\+D\n){n\+A\n + l){n\, (46) 



l<:c<iV + l 



where 

/ \ f 2^ {x < n), 

' \ X - \ (x > n). ^ ' 

By using the relation DA = qAD, one can check that this matrix satisfies the following relation, 
a "degenerated version" of the commutation relation (PO]) : 

^jW(„(JV,„) ^ ^iN,n)^ _ (^(w,„) ^ a(^-»))M(;y-i) = 0. (48) 
From this relation, the L-fold tensor product of a'^^") satisfies 

l<i<L-l ' KjXL-1 
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Noting that the summations of local Markov matrices are total Markov matrices, we get 



where = 



AN, 



v)). Thus we find the submatrix V's,s\{s„} of vj/C^^")^ i.e. 



V's,6\{s„} : V^\{s„} ^ {ji---jL\-^^,s\{s„}\ki---kL) = {{w\a 



■ a I. \v) 



with a 



{N,n) 
jk 



(j|a(^'")|fc), satisfies the conjugation relation 

^sV'5,s\{s„} = i^!i,5\{s„}Ms\{s^}. 



(50) 



(51) 



(52) 



Each element (|5f p of the conjugation matrix ?A5^s\{s„} becomes or a power of q with its exponent 
corresponding to how many times we need to exchange DA — > AD for reordering the matrix product 
to A ■■■ AD ■■■ D: 



Sn—Sn — 1 



{ji- ■ ■ jL\-4^<i,5\{s,,}\kl- ■ -kL) 



exp (ing ^ < = n}^ {Xn{ji) = fc< (Vi)), 



(53) 



(otherwise). 
For example, for the sectors s = {1, 3, 5} O (1, 2, 2, 1) and s \ {52} = {1, 5} O (1, 4, 1) with L = 6, 
(123423|V's,;,\{s2} 1122322) = = {{w\DADA\v)) 



q{{w\DAAD\v)) = ^^^(uilADADlt;)) = {{w\AAD D\v)) 



(54) 



The generalized matrix (product) ansatz, i.e. the procedure (P5|) - ([5^ . was introduced in [3] for 
the periodic boundary condition. There the right-hand side of (|48)) is replaced as 



with an auxiliary matrix a*-^'"\ and the conjugation matrix is constructed by taking the trace 



{ii ■ ■ •iL|^s,5\{.s„}l^i • • • fcL> = Tr 



{N,n) 



(56) 



(We set a'^^'"^ = in our case.) The families of representations for the hat relation (1551) found in [HIS] 
are more complicated than our case. For example, for (N, n) (3, 1), our representation (I46p is 



2 3 



1 /D 0\ 



,(3,1) 



A 
10 

V ly 



(reflective). 



(57) 



On the other hand, the representation found in 4 , 5 is 



where 



,(3,1) 



f l(g>l (5(8)1 l(g)S\ 

A(^A 

e(E)A 1(E) A 

yi(Xie (5(g)e l(g)iy 



/ ci 

C2 



V 



/ 

ci 

C2 

V 



(periodic). 



(58) 



(59) 
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and a^^'") needs to be chosen as diag(l, q, q, q)a'^^'-^\ Note that the number of tensor products of each 
element in the solution for the periodic case increases as N increases. On the other hand, the solution 
(|46l) does not contain a tensor product. 

The hat relation (jSSp is independent from boundary conditions. However, whether a representation 
for the algebra defined by the hat relation is practical (i.e. whether a representation allows us to 
construct a nontrivial conjugation matrix) depends on boundary conditions. For example, the matrix 
V's,s\{s„} defined by with the representation ([571) is 0. 

In the same way, we can construct the restricted version for (p psp . starting with the matrix 



a 



(N,n) = J2i<:^<N+i l\xnix)){x\+D\n){n+l\. The matrix defined by (fci ••• •• -jl) 

-^s\{s„}Vs\{s„}s 



This matrix has indeed trivial elements 

'Pb\{s„},s = ■ ■■XnUL)){jl ■ • -JlI 



satisfies the conjugation relation 'y3s\{s„}s-^^s 



(60) 



where the summation runs over all the configuration in the sector s. For ket vectors, this matrix 
"identifies" nth and {n + l)st class particles as a same class [6^. The matrix V's,s\{s„} (ED) sends 
any basis vector \ki ■ ■ ■ k^) e Vg\{s„} to vectors in such that each particle ki keeps its position 
i, (s„+i — Sn) particles of nth class are changed to (n + l)st-class particles (the rest of (s„ — Sn-i) 
particles of nth-class are unchanged), and i^th-class particles {v > n) are changed to {v + l)st-class 
particles. Such vectors are restored to |fci • • • ki) by the identification matrix '^s\{s^},s'- 



^s\{s^},^i^s,s\{s„}\kl ■ ■ ■ kL) ^V'At^.J.s X! ((^lOjifci ■■■0-jLkL\v))\il---iL) 



X7i(ji) = fci (Vi) 



Y iHa'Jiki---ajr^k^\v))\ki---kL) =: C\ki ■ ■ ■ kL) 



Actually the constant C is independent of the configuration ki ■ ■ ■ k^ and one can show 



C = 



E 



{{w\Ui---Us. 



■. + l-s„ 



Thus we have 



<PS\{S„},5V'S,S\{S„} 



l]![Sn-|-l 



(61) 



(62) 



(63) 



where Id is the identity matrix. The injectivity of V'5,5\{s„} follows from this relation, and we have 
the inclusion relation 

SpecMs D Speci\4\{s„}- (64) 
Now we turn to the construction of the conjugation matrix between Markov matrices of N- and 
iV'-species sectors {N — N' — u > 0). Let us consider the local relation, which we also call hat relation. 



0. 



(65) 

We have already known a family of solutions to this relation: 

A' = a(^'"i)*...*a(^'+i'"") (1 < < A^-^+ 1), (66) 

*9«i ® ^ife| . matrix- valued matrices 



where the symbol * denotes the product Q ^ R 
Q — {Qij}ij and R = {Rij}ij- For example, 

/ 1 \ 

D 

A 

\ 1 J 



X 




( t®D 

D®A 

A® A 

V 1(8)1 j 



(67) 
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is a solution to mII{{X X) = {X ® X)mIII. 

Bookending each element of A"®^ between {{w\®" and \v))^^, we obtain 

* = ((wl^^A-^^lw))®" = • • • ^(^'+1.""), (68) 

which satisfies 

AfW* = *A/(^'). (69) 

The matrix ^f* sends a vector of an A^'-species sector to an A^-species sector, via (iV' + l)-species — )■ 
(iV' + 2)-species {N — l)-species. Each index ng specifies which class of particles splits in 

sending an (N — ^)-species vector to an {N — £ + l)-species sector. 

For the A^-species sector s = {si < ■ ■ ■ < s^} and the A^'-species sector t = s \ {s^-^, . . . ,Si,^^} 
{u = N — N'), we have the conjugation matrix 

i'si = V's,s\{s„i}V'0\{s„J,s\{s„^.s„2} ■ ■ •'0tu{s„„},t (70) 

satisfying 

M,^M = ^,tMt . (71) 

This is a submatrix of ^ (j68p with the choice ni ~ vi — =^{z\z < £,i>z < i^i}, but indeed independent 
of the choice. It is enough to show the simplest case u — 2: 

"comniutativity" : V'0,s\{'5f.}V'A{v}^A{V:S^} = V's,s\{s„}V'0\{s„},5\{V:S.4- C^^) 

We suppose ^jl < v, and use the explicit expression ((5^ . The choices of s for the left- and right-hand 
sides of ([7^ are [111,112) — — 1) and respectively. By the expression ([55]) . both sides are 

calculated as 

{jl ■ ■ ■ jL\i^^^,\{s^}^,\{s^}^,\{s^^s,^}\h ■ ■ ■ kL) = (73) 



= h (Vi) 



exp(ln(j( ^ #{«'<«|>=m}+ Y1 #{*'<«IXM(ji') = J^-l}^ 

(otherwise) , 

{ii---iL\'4^<,,^\{s^}i^<i\{s,,},<,\{s,,,s,,}\ki---kL) = (74) 



fci (Vi) 



exp(lnq( < ^ v} + H^{i' < i\Xu{3^') ^ ^ 

(otherwise), 
which are equal. 

5 Stationary state 

The zero-species sector consists only of the configuration 1 • • • 1. Since M'f^^ = 0, the hat relation 

(|65| with A^' = becomes m\^1{X ® X) ^ Q. A solution to this hat relation is given as X = 
^(AT.ni) ^ ^(Af-i,n2) ^ . . . ^ a^^'^\ and thus we obtain a stationary state in the matrix product form 
[T1[71[TT]: the probability P{ji ■ ■ ■ jl) of finding a configuration ji ■ ■ ■ Jl can be expressed as 



Pin ■ --jl) = ^{{w\^''X,, ■ • •x.jT.r^, (75) 



where 
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satisfying XaXp — qXpXa {a < 13) 0. This implies that the system satisfies the detailed-balance 
condition as we commented in section [21 

The stationary states in the periodic boundary condition can also be written in the matrix product 
form [H [19] . However the representation for the matrices {Xi 's) in our case is much simpler than that 
of the periodic-boundary case. For example, for = 3 with the choice n2 = ri3 = 1, we have 



a(3.i)^a(2,i)^„(i,i) 



The stationary state |Po) 
matrices and |1 • • • 1) as 



( X,\ 

X2 

Xs 

V ^4 J 

/ 1(8)1 



/ D®D®D \ 

1(g) A(g) D 
\ t®t®A J 

{1 + 5) + 5 <^l(Si A + \(^ 5 ® {e +1) \ 

A®A®{t + 5) 
e®A®{l + 5) + l®A 
{t + 5) + 5®e®A+ti 



(reflective) 



A 



V 1(8> 
of each sector s = {si < 



(77) 



(periodic). (78) 



l(8(e-t-l) ) 
< sjv} is given by a product of conjugation 



(79) 



with {vi, . . . ,vm} = {1,...,-/V}. In other words, the stationary state in sector s is constructed 



by transferring the zero-species vector to via sectors {siy„} — >■ {si^jvi ■^I'jv-i} 



Note that the stationary state (|7^ is indeed independent of the choice of intermediate sectors (i.e. 
independent of the choice of vi s.) Noting equation we find Z of the sector s = {si < • • • < sjv} 
in the general case is given by the g-multinomial 

= (1 • • • l|</'0{s„} • • • '^s\{s^],s i^B,s\{s^} ■ ■ ■ '0{s„}0|l • • • 1) 

[L]\ 



(80) 



[si - sq\\[s2 - Si]! • • • [sN+1 - sjv]! 



(so = 0, sjv+i = L), 



which is the dimension of the sector s for q — 1. 

For the TASEP case q — 0, a. unique configuration can be realized in the stationary state since 
XaXis = (a < /3). That is, the stationary state in each sector is an absorbing state, where all the 
particles stay in the descending order ji ■ ■ ■ Jl [ji > ji+i for Vi). On the other hand, for the SSEP 
case (? = 1, all the possible configurations are realized with an equal probability. 

Now we turn to the relaxation to the stationary state, where the relaxation time t is characterized 
by the largest non-zero eigenvalues as r = — Rei?~^. (The largest eigenvalue is indeed 0, which 
corresponds to the stationary state.) We first consider the simplest sector {L — 1} o {L ~ 1, 1), i.e. 
L — 1 partiles and 1 vacancy. The spectrum of M{i_i} is given by Spec(A/{i_i}) {0} U { — (1 -f g) -I- 
2y/gcos -j^\k In particular the largest non-zero eigenvalue is 



E 



-(1 + 9) 



27gcos-. 



(81) 



The spectrum of the general one- species sector {si} O (si,L — si) (1 < Si < L) contains that of the 
one- vacancy sector {L — 1}. Furthermore we expect that the largest non-zero eigenvalue of the sector 



The stationary state for the multis-pecies ASEP with the general hopping rule a/3 — > /3« (rate F^^) has also the 



matrix product form with the algebra XaXp 



-XpXa, see [ij for N = 3. We need a 



-^^i^^il^-fold tensor product 



in the representation for this algebra, and the decomposition structure II76I I no longer exists. 

^ This can be easily derived by the Bethe ansatz, and eigenvectors (except the stationary state) are given by 



\Pk){L-l} = I2l< 



\1- ■ -1 2 1- ■ -1) with A : 



1, 



!)• 
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{si} is equal to that of {L — 1} ([5T|) . which can be checked in small systems. This observation implies 
that the relaxation time behaves as 

r ^ (1 - ^)-^ {q < 1), —{q = 1), (82) 

TT 

as L — )• c». The "fist excited state" \Pi)_^g^-j, i.e. the eigenvector corresponding to the largest non-zero 
eigenvalue of the sector {spf} can be written as |-Pi){sjv} = '0{sjv}{i-i}l-Pi){i-i} where '0{sjv}{i-i} 
constructed as ((37|) . Recall that the spectrum of the general multi-species sector s = {si < ••• < sn} 
also contains that of the sector {sjv}, see ((64)) . Again we expect the largest non-zero eigenvalue of the 
sector s is identical to (|5T|). which can be checked for sectors with small dimensions. This implies the 
same behavior of the relaxation time as for the one-species case. The corresponding eigenvector 
of the sector s also has the form |Pi)s = ■*/'s.5\{s„j } ' ' ' V'{sn-i,sjv},{sjv}I^i){sjv}- I'^ contrast to our case, 
the relaxation time of the multi-species ASEP behaves as 

r ^ Li (9 < 1), (q = 1) (83) 

in the periodic boundary condition [SI [T3] . Comparing the behaviors (|82p and ([55)) , we notice that 
the boundary condition plays an important role. 

6 Summary 

We investigated a multi-species generalization of the ASEP with reflective boundaries. We found the 
symmetry of the Markov matrix can be interpreted as a matrix product form, constructing conjugation 
matrices (|57)) and (|55)) which intertwine Markov matrices of different sectors. We showed that the 
conjugation relations follow from the local commutation relations (HO))- We also considered relations 
between Markov matrices of different values of N. We constructed a conjugation matrix connecting 
dynamics of a simpler system and a more complex system by using solutions to the hat relation ()55p . 
We saw that the stationary state can be written in a product of conjugation matrices, and the first 
excited state can be obtained by multiplying that of one-species sector by a product of conjugation 
matrices. These properties are also true in the periodic boundary condition. (However the behaviors of 
the relaxation time to the stationary states are different in these two boundary conditions in general.) 

It is remarkable that there exist several solutions (representations) a^-^'"^ to the hat relation (|55)) 
which are suitable either the reflective boundary condition or the periodic one [H [S] , and the choice 
of the hat matrix a^^'"'^ have to be changed according to the boundary conditions. At present we 
do not have a systematic way to flnd appropriate representations as well as appropriate hat matrices. 
Another interesting study will be applying (or generalizing) our method to the system with injection 
and extraction of particles [71 |Tl] . 
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